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ABSTRACT 


Mass is one of the most fundamental parameters characterizing the dynamics of a coronal mass ejection (CME). It has 
been found that CME apparent mass measured from the brightness enhancement in coronagraphs increases during its 
evolution in the corona. However, the physics behind it is not clear. Does the apparent mass gain come from the 
outflow from the dimming regions in the low corona, or from the pileup of the solar wind plasma around the CME? 
Here we analyze the mass evolution of six CME events. Based on the coronagraph observations from the Solar 
Terrestrial Relations Observatory, we find that their masses increased by a factor of 1.3—1.7 from 7 to 15 Rs, where 
the occulting effect is negligible. We then adopt the “snow-plow” model to calculate the mass contribution of the 
piled-up solar wind. The result gives evidence that the solar wind pileup probably makes a non-negligible 
contribution to the mass increase. In the height range from about 7 to 15 Rs, the ratio of the modeled to the measured 
mass increase is roughly larger than 0.55 though the ratios are believed to be overestimated. It is not clear yet whether 
the solar wind pileup is a major contributor to the final mass derived from coronagraph observations, but it does play 


an increasingly important role in the mass increase as a CME moves further away from the Sun. 
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1. INTRODUCTION 


Mass is a major parameter characterizing the dynamics of a 
coronal mass ejection (CME). However, due to the limitation of 
observations, we do not know very much about the CME mass, 
and its evolution from the corona to interplanetary space. CME 
masses are mostly derived from images of white-light 
coronagraphs (e.g., Vourlidas et al. 2000), as a CME can 
cause an enhanced brightness in white-light images. However, 
the mass derived from white-light coronagraph images is just 

n “apparent” mass of a CME. It has two components: one is 
the mass of plasma really ejected from the corona, the other is 
the piled-up solar wind mass due to the compression of the 
solar wind plasma surrounding the CME. It is difficult to 
separate the two components from imaging data. Therefore, 
what we often calculate is the apparent mass. A few studies 
have shown that the apparent mass of a CME increases with 
time (e.g., Colaninno & Vourlidas 2009; Bein et al. 2013). 

The statistical work by Wang et al. (2011) has shown that the 
CME brightness is positively correlated with the CME speed. A 
natural explanation is that a faster CME causes higher 
compression of ambient solar wind and therefore a brighter 
signature in coronagraph images. Numerical simulations also 
revealed that the mass of the high-density shell ahead of an 
erupted flux rope increases as it propagates away from the Sun 
but the mass of the dark cavity, where the flux rope is located, 
remains almost unchanged (Lugaz et al. 2005; Riley 
et al. 2008). All these studies imply that the piled-up solar 
wind mass may not be negligible, and may contribute 
significantly to the aforementioned phenomenon—apparent 
mass gain of a CME during its propagation. 

On the other hand, there is also evidence suggesting that the 
CME mass gain can be attributed to the continuous outflow from 
the low corona (e.g., Howard & Vourlidas 2005; Tian et al. 
2012; Bein et al. 2013). Bein et al. (2013) investigated the CME 


mass evolution for a set of 25 events, and found that the apparent 
mass increased by about 2%-6% per solar radius based on the 
COR2 observations taken by Solar Terrestrial Relations 
Observatory (STEREO) (Howard et al. 2008; Kaiser et al. 
2008). They further found that the mass centroids of these CMEs 
moved toward the Sun relative to their self-similar expansion in 
most events. The authors suggested that the reason might be that 
the rear part of the CME with higher mass might propagate 
slower than the front part of the CME. They found that the 
higher mass at the rear was supplied by the outflow from the 
lower corona. However, in the analysis of Bein et al. (2013), the 
fact that materials may fall back to the Sun, which can also cause 
a backward motion of the CME centroid, was not considered. 
Thus, the question how much pileup and outflow contribute to 
the mass increase still remains unsolved. 

In this paper, we are going to answer the question why the 
CME apparent mass increases in a way different from that by 
Bein et al. (2013). We try to quantitatively figure out at which 
height and how significantly the piled-up solar wind mass 
contributes to the CME apparent mass. In Section 2, we briefly 
describe the selection of events, their observations, and their 
three-dimensional (3D) geometric and kinematic properties. 
How we derive the CME mass and the mass evolution of 
selected events is presented in Section 3. In Section 4, we 
analyze the possible source of the CME mass with the help of a 
solar wind pileup model and MHD numerical simulations. The 
final section consists of the discussions and conclusions. 


2. EVENTS 
2.1. Event Selection 


Six CME events (see Table 1) were selected from the halo 
CME list (http://space.ustc.edu.cn/dreams/fhcemes/) com- 
piled by Shen et al. (2013a). The criteria to choose these 
events are: (1) they appeared halo in the Large Angle 
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Table 1 
The Geometric and Kinematic Parameters of Six CME Events 


No. Date oe) @(°) 
El 2010 Feb 12 —2.5 Ta 
E2 2010 Apr 03 4.6 —25.1 
E3 2011 Apr 17 -177.1 8.9 
E4 2011 Jun 21 —20.8 9.5 
E5 2011 Aug 03 15.4 16.4 
E6 2012 Feb 23 67.2 24.4 


We (°) W. ©) Q (sr) Vavg (km s=») 
68 47 1.0 752 + 31 
TH 50 1.2 824 + 29 
72 52 1.1 1312 + 51 
70 70 1S 1006 + 38 
68 41 0.8 1294 + 56 

139 82 3.5 445 + 12 


Note. ¢ and @ are the longitude and latitude of the 3D CME propagation direction, respectively. Positive/negative values of ¢ mean westward/eastward propagation 
with respect to the Earth. Wp = 2(a + ô) and W, = 26 are the face-on and edge-on angular width, respectively. Q = Wp W, = 46(a + 6) approximates the solid angle 
of a CME. vavg is the average speed of the CME from the linear fit to the observed height—time profile. 


Spectrometric Coronagraph (LASCO) on board the Solar and 
Heliospheric Observatory (SOHO) spacecraft. They were also 
observed by COR1 and COR2 coronagraphs on board STEREO 
and in most cases as almost limb events. Limb CMEs are less 
occulted by the coronagraph occulter, thus are favorable for 
mass estimate. (2) There was no contamination from other 
CMEs in the region where the investigated CME was located. 
(3) There was no contamination of very bright Ha line 
emission from prominence material inside the investigated 
CME in the coronagraph images. Mass estimate of such 
contaminated CMEs is often misleading if Thomson scattering 
theory is applied (Mierla et al. 2011; Carley et al. 2012). The 
six CMEs cover a broad range of leading-edge speed from 450 
to 1300kms_' as indicated in Table 1. 

The twin spacecraft STEREO A+B and SOHO monitor the 
solar corona from three different angles of view at about 1 AU. 
COR1 and COR2 on board STEREO observe the corona from 
1.4 to 4 Rs and 2.5 to 15 Rs, respectively (Howard et al. 2008). 
LASCO C2 and C3 cameras on board SOHO have a field of 
view (FOV) of 2—6 Rs, and 3.7—-32 Rs, respectively (Brueckner 
et al. 1995). 


2.2. Deprojected CME Geometry and Kinematics 


The projection of a CME onto an image may significantly 
distort its geometric and kinematic parameters. For a better 
estimate of the mass of a CME and its evolution, we have tried 
to derive the deprojected 3D geometry and kinematics. Various 
models/methods have been developed to retrieve realistic 
parameters based on multi-point imaging data, e.g., the forward 
modeling (e.g., Thernisien et al. 2009; Wood et al. 2009), 
triangulation method (e.g., Liewer et al. 2009; Temmer 
et al. 2009), polarization ratio method (e.g., Moran 
et al. 2010), inverse method (Antunes et al. 2009), geometric 
localization (e.g., de Koning et al. 2009; Byrne et al. 2010; 
Feng et al. 2012). The comparison of various models could be 
found in, e.g., Mierla et al. (2010), Thernisien et al. (2011) and 
Feng et al. (2013a). 

Here we use the forward modeling developed by Thernisien 
et al. (2009) to derive the 3D geometric and kinematic 
parameters of the six CME events. In this method, the 
graduated cylindrical shell (GCS) model is applied, which 
consists of a tubular section forming the main body of the 
structure attached to two cones that correspond to the “legs” of 
the CME. This model has six free parameters: the propagation 
longitude and latitude (¢ and 0), aspect ratio («), tilt angle(y), 
half angle between two legs (a), and heliocentric height of the 
CME leading edge (A). From the parameters a and k, we 
further derived the edge-on angular width (26) and face-on 


angular width (2(a@ + 6)), where 6 = arcsin(«) is the half 
angle of the cone, i.e., the CME leg. Note that all the heights 
mentioned in this work are the radial distance measured from 
the solar center. These six parameters were tuned to obtain a 
best match between the modeled flux rope and the observed 
CME from different viewpoints. More details of this model and 
its application to CMEs could be found in, e.g., Thernisien 
et al. (2009, 2011) and Gui et al. (2011). For the time sequence 
of images of each CME, we obtained a set of parameters at 
each time. For simplicity, the free parameters are set to be 
constants by trial and error, except the height h which is a 
function of time. 

As a demonstration, Figure 1 shows the results for the event 
on 2010 April 3. The upper panels are the images taken by 
COR B, C2 and COR A, respectively, in which a pre-event 
image has been subtracted to remove the background and make 
the CME structure more pronounced. To get a picture of the 
CME as complete as possible, we incorporated CORI and 
COR2 observations. The aforementioned procedure has been 
applied to all the other five CME events as well. The lower 
panels are the same as the upper panels except that the green 
meshes are superimposed to show the best-fit flux rope. It has 
to be pointed out that for the GCS reconstruction, we did not 
include in the fitting process the shock—-sheath region when 
visible. By repeatedly applying the GCS model to the time 
sequence of the images of the CME, the evolution of the height 
of the leading edge of this CME can be obtained and is shown 
in Figure 2. The uncertainty in the height is indicated by the 
error bars. It has two sources; one is the ~0.3 Rs uncertainty in 
the definition of the CME leading front in the white-light 
images (e.g., Bein et al. 2013), and the other comes from the 
forward fitting process. We assumed that the uncertainty is 
about 5% of the height at each time instance. At the height of 
20 Rs, it produces an uncertainty of 1 Rs which is close to the 
maximal uncertainty in the sensitivity analysis in Thernisien 
et al. (2009). The height-time data points were further fitted by 
a linear function to estimate the average speed of the CME. 

Table | lists the geometric and kinematic parameters of our 
six events, i.e., longitude and latitude of the propagation 
direction (¢ and 0), face-on and edge-on angular width (Wp and 
W.), CME solid angle (Q), and average speed (Vayg). The 
longitude and the latitude are in a heliographic coordinate 
system with Earth sitting at zero longitude. Positive /negative 
longitude represents a CME propagating westwards /eastwards. 
The solid angle spanned by a CME is approximated by the 
multiplication of the face-on and edge-on angular width. The 
leading-edge speeds of the CME events vary from about 450 to 
1300 km s~, showing good representativeness. 


THE ASTROPHYSICAL JOURNAL, 812:70 (12pp), 2015 October 10 


FENG ET AL. 


Figure 1. Upper panels: the CME observed by STEREO COR1 and COR2 and by SOHO/LASCO C2 on 2010 April 3. Note that COR1 and COR2 images are 
incorporated to have a comprehensive overview of the CME morphology. Lower panels: the same as the upper panel except that the modeled flux rope indicated by 


green meshes are superposed. 
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Figure 2. Deprojected height as a function of time for the CME event on 2010 

April 3. The asterisks are the measurements and the black solid line is a linear 


fit. The average speed derived from the linear fit is marked in the upper part. 
The error bars indicate the 5% uncertainty of the deprojected heights. 


3. EVOLUTION OF CME APPARENT MASS 
3.1. Mass Estimate Method 


Mass calculations based on white-light coronagraph images 
are usually based on the Thomson scattering theory (Min- 
naert 1930; van de Hulst 1950; Billings 1966; Howard & 
Tappin 2009). The major non-Thomson scattered contributions 
like stray light and dust scatter are removed by subtracting a 
pre-event background image. The remaining image brightness 


then should show the CME density distribution produced by 
Thomson scattering alone. For the total brightness in a given 
image pixel, the details of this scattering process can be 
described by 


T 


Oe 
Ti F (r, x, u) Ne(r) de, 
(1 — wC(r) + uD(r) 


F(r, x, u)=2 1-473 
(1 — w)A(r) + uB(r) . 5 
1—u/3 sin’ x (1) 


where bo is the physical mean solar brightness (MSB) that is 
used as a unit of the pixel value in calibrated CORI and COR2 
observations, F(r, x, u) is a function of three variables r, x, and 
u. r is the distance of the scattering location to the solar center, 
x is the scattering angle between the radial vector from the 
solar center to the scattering electron and a position vector from 
the electron toward the observer. u accounts for the solar disk 
limb darkening, and we used the value of 0.56 for u in the 
visible white-light spectral range. A, B, C and D are known 
functions of r. a, is the differential Thomson scattering cross 
section, and has a value of 7.95 x 10-76 cm? sr™!. 

In Equation (1), particular attention should be paid to the 
scattering angle y. In this work, the scattering angles were 
computed under the assumption that a CME is concentrated in 
the plane defined by the CME propagation direction and the 
rotational axis of the Sun as in Colaninno & Vourlidas (2009) 
and Carley et al. (2012). And the CME propagation direction 
was derived in Section 2.2 using GCS forward modeling 
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Figure 3. Mass images within the CME sector regions. The plus sign and a white circle indicate the solar center and solar disk, respectively. The mass color bar is 


shown on the right side. 


method. Under such all-in-propagation-plane assumption of a 
CME, we can calculate the column density J òs N, dl for any 
LOS with Equation (1). Assuming a composition of 90% H and 
10% He (Vourlidas et al. 2000), the mass in each pixel of a 
coronagraph image can be derived by multiplying the column 
density with the pixel area measured at a distance of about 
1 AU from the observer. 

In our practical procedure, we first defined a sector that 
contains the whole CME. In Figure 3, the mass in each pixel 
within the sectors in COR A and COR B is presented for the 
CME occurred on 2010 April 3. And we integrated the mass 
over the entire sector and obtained two mass values from COR 
A and COR B data, respectively. It should be noted that as we 
subtracted a pre-event background from the image of interest, 
the brightness in the non-CME region is supposed to be 
roughly at noise level. Thus, although the sector is larger than 
the region that the CME actually occupies, the calculated 
values of the apparent mass should be acceptable. Due to the 
different perspectives of COR A and B, the two values are 
slightly different. In Figure 4 we show the mass evolution with 
time. The mass evolution derived from COR A observations is 
delineated by the red solid line, and from COR B observations 
by the red dashed line. The average of the mass derived from 
COR A and B observations is delineated by the green solid 
line. We found that the difference between the mass derived 
from COR A and B was about 10 + 4%. Some other curves 
and associated error bars in Figure 4 will be explained in 
Section 3.2. The evolution of the averaged mass derived from 
COR A and B data for all six events is shown in Figure 5. The 
mass uncertainties will also be described in Section 3.2. 


3.2. Mass Uncertainty Analyses 


Vourlidas et al. (2010) has shown that the instrumentally 
photometric uncertainty is quite small, and the subtraction of a 
pre-event image introduces insignificant errors, if the pre-event 
image is selected carefully. Similar to Feng et al. (2015), we 
mainly considered two sources of error in the mass estimate. 
One is the LOS distribution of the CME mass in Equation (1), 
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Figure 4. Mass evolution with time for the CME on 2010 April 3. The thin 
black solid lines are the mass derived from COR A total brightness images with 
different selections of CME sectors and varing propagational longitudes within 
$ + 5°. A total of 15 independent mass measurements were made for the error 
estimate. The mass along the thin black dashed lines are derived in a similar 
way but from COR B observations. The mean mass evolution derived from 
COR A and COR B are marked by solid and dashed lines in red, respectively. 
Fifteen measurements from COR A and 15 measurements from COR B yield 
15 averaged mass over A and B. The green line and associated error bars are 
the mean and 30 values of the 15 averaged mass measurements. 


especially the longitudinal propagation angle œ. The other 
source is our somewhat arbitrary selection of the CME sector. 
To estimate the uncertainties from these two sources, we varied 
ġo within ¢ + 5° at seven values and repeated the CME sector 
selection process independently for eight times. A total of 15 
measurements were made from COR A and B observations, 
respectively. They are shown by black solid and black dashed 
lines in Figure 4. Correspondingly, the red solid and dashed 
lines are the averaged values over 15 measurements. For our 
analyses of the mass evolution of a CME event, we adopted the 
green curve. Fifteen measurements from COR A and 15 
measurements from COR B yield 15 averaged mass over A and 
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Figure 5. Mass evolution of all six events. The black error bars are the 3c of all the 15 mass measurements with different sector regions and different propagation 
angles. The red error bars are the 30 of eight mass measurements only with different sector regions. For most cases, major error comes from the selection of the CME 
sector. While for the CME on 2012 February 23, the major error comes from the uncertainty of the propagation angle. 


THE ASTROPHYSICAL JOURNAL, 812:70 (12pp), 2015 October 10 


ae ee Ta 

5 L A 
S 0.8 `A J 
m [ \\ 

j f \\ 

x< E \\ 

a 0.6f \ 4 
ie ‘ 

SS | \ al 
= 0.4 F Ñ 

3 \ 

5 L \ 

>< 0.2 \ a 
a E R=10 RÀ 
> L R=4 Rg 
OM in 6s a ee ge ee 

QO 50 100 150 
x (degree) 


Figure 6. Dependence of the function F(r, x, u) on the scattering angle x at 
three impact distances (R = 2, 4, 10 Rs). For COR observations, the impact 
distance is the projected distance of r from the scattering location to the solar 
center onto the POS. The brightness is normalized by its maximal value at the 
POS (F(r, x = 90°, u)). The vertical dotted line marks the position of x = 90°. 


B. The green line and associated error bars are the mean and 30 
values of the 15 averaged mass measurements. 

Figure 4 presents an example of the mass evolution and error 
analyses. In Figure 5, the mass evolution and associated error 
bars of all six events are shown. The black error bars are the 3a 
of all 15 mass measurements with different sectors and 
propagation longitudes, while the red error bars are the 3c of 
the mass measurements only with different sector regions. In 
most cases the major uncertainty comes from the selection of 
the CME sector. However, for the CME observed on 2012 
February 23, the major uncertainty originates from the error of 
the propagation longitude. 

To explain the results of uncertainty above, we include in 
Figure 6 the dependence of the function F(r, x, u) on the 
scattering angle x at three impact distances of 2, 4, and 10 Rs. 
The ordinate is F(r, Xx, u) normalized by its maximum at the 
POS x = 90°. For COR observations, the impact distance is the 
projected distance of r from the scattering location to the solar 
center onto the POS. As revealed by Figure 6, when x deviates 
from 90° within about 30°, F(r, x, u) only changes slightly and 
forms a plateau around x = 90°. This phenomenon has also 
been found by Howard & DeForest (2012) where they termed it 
“Thomson plateau” for the heliospheric HI observations and 
Colaninno & Vourlidas (2009). Therefore, a small change in 
propagation longitude ¢ + 5°, which causes a small change of 
x + 5°, does not effect F(r, x, u) much within the 30°- 
deviation range. Subsequently, the column density and mass 
for a given pixel derived from Equation (1) does not change 
much when y deviates from 90° within about 30°. This is the 
case for the first five CME events. Therefore, the uncertainty of 
the mass for these events mainly comes from the selection of 
the CME sector rather than the propagation longitude. 

In the case of the CME occurred on 2012 February 23, the 
propagation direction deviated from the POS as seen by 
STEREO A by about 50°, and from the POS as seen by 
STEREO B by about 90°. As the CME appeared in the FOV of 
COR B as a full halo, the corresponding mass calculation 
subjects to a large error. Therefore, we only used COR A data 
for mass analyses. When we derived the mass from COR A 
data, a small change in y yielded a large change in mass. It is 
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because the slope of F(r, x, u) at x ~ 140° is quite steep in 
Figure 6, a small change of x + 5° makes a large change in F(r, 
x, u), hence in column density and mass. That is why the major 
uncertainty of this event comes from the uncertainty of the 
propagation longitude. 

Note that our all-in-propagation-plane assumption yields an 
infinitesimal angular width for a CME. Mass can be under- 
estimated under such an assumption according to Vourlidas 
et al. (2000). The authors derived the ratio of the mass estimate 
under the infinitesimal-angular-width assumption to the mass 
under a finite-angular-width assumption as a function of the 
angular width of a CME. As the angular width usually keeps 
nearly constant in time in the COR FOV, the under-estimate 
ratio should also keep constant. Therefore, the mass change rate 
dm/dt is not effected by the infinitesimal-angular-width 
assumption. 


3.3. Treatment of the Occulting Effect 


Figure 5 clearly shows that the mass evolution has an 
increasing trend. Bein et al. (2013) have demonstrated that one 
of the major causes of such an increase is the entering of the 
CME ejecta into the COR FOV from behind the occulter, 
which we call occulting effect. They then showed that the 
apparent mass evolution could be fitted by the following 
function 


3 
m(h) = mf — (+=) | + Am(h — hocc), (2) 


where the parameter hocc is the height of the outer edge of an 
effective occulter, mo is a modeled initial apparent mass of a 
CME, and Am is an additional mass increment per unit height. 
The first term on the right-hand side of Equation (2) represents 
the increase of visible mass from a homogeneous cone-shaped 
and self-similarly expanding CME. It describes how the visible 
mass of a CME with a constant mass mọ evolves, when 
entering the FOV of COR due to its expansion. The increase of 
the visible mass above the occulter is merely a geometric effect. 
The second term is an empirical correction to this simplified 
model evolution, as the CME mass could increase due to 
physical reasons, e.g., mass outflow from a dimming region, or 
solar wind pileup. In the model of Bein et al. (2013), the actual 
physical mass increase of a CME can be approximated by 
m(h) = mo + Am(h — hocc). 

In Figure 7, the black discrete data points are the calculated 
mass as a function of deprojected height with the error bars in 
horizontal and vertical directions. Taking these error bars into 
account, we fitted Equation (2) to the mass-height measure- 
ments of each CME event. The best-fit parameters for each 
CME event are written in the bottom right, and for mp and Am, 
the mean +ø are in logarithmic scale. In each panel the fitting 
results are shown by the black solid line, and the black dashed 
line is the mass variation mo + Am(h — hoc) with the 
occulting effect removed. The vertical line marks the height 
from which the occulting effect starts to be negligible. It is 
defined as the height from which the measured mass can reach 
above 97% of the dashed line. It is found that the occulting 
effect could be ignored once the CME reached about 7-8 Rs. 
We denote this height as h4, the corresponding mass at h; is mı. 
mı and h; for all six CME events are listed in Table 2. The 
other two mass quantities that are also included in Table 2 are 
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Figure 7. Mass profiles as a function of the deprojected height. The best fit result of Equation (2) is indicated by a black solid line in each panel. The mass evolution 
after removing the occulting effect is shown by the black dashed line. The best fit parameters are marked in the bottom right in each panel. For mg and Am, the mean 
+o is written in logarithmic scale. The vertical dotted line marks the height A; where the occulting effect starts to be negligible. The mass evolution due to the solar 
wind pileup above h, is shown by a red dashed line in each panel. The uncertainties indicated by the shadow region in orange are derived under the assumption that 
Ny and v,,, vary in the range from 80% to 120% of the simulated values. 


mo below hocc, and the final apparent mass mz before a CME negligible. We did not use this simplified model for further 
leaves COR FOV at hp. analyses of mass evolution. The measured mass evolution may 

Note that our purpose of adopting this model is only to deviate from the linear increase assumption in Equation (2). 
obtain the height where the occulting effect starts to be Actually, in Figure 7, for most events, when the CME starts to 
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Table 2 
Derived Parameters of CME Mass at Different Heights, Solar Wind Parameters, and its Mass Contribution 
No. mo hoce mı hı m hz New Viw my f 
logio(g) (Rs) logio(g) (Rs) logio(g) (Rs) (cm ~*) (km s~’) logio(g) 
El 15.17 2.83 15.27 7.74 15.37 14.2 869 263 15.09 £2.93 
E2 15.38 2.34 15.48 7.16 15.60 14.3 657 376 15.01 <1.06 
E3 15.35 1.82 15.63 7.01 15.86 15.7 824 243 15.32 <0.70 
E4 15.45 2.69 15.62 7.97 15.78 15.8 696 355 15.20 <0.85 
E5 15.42 2.62 15.68 8.46 15.85 15.6 878 226 15.15 <0.62 
E6 15.65 2.79 15.83 7.51 16.01 15.3 610 338 15.28 <0.55 


Note. mo is the modeled initial apparent mass of a CME below the effective occulter height hocc. mı is the apparent mass at the height h, where the occulting effect 
begins to be negligible. m is the final apparent mass before a CME leaves COR FOV, and h, is the corresponding height. nsw and v,, are the ambient solar wind 


density and speed at h = 20 Rs. m, is the solar wind mass piled up around the CME from h; to ho. f, = 


the apparent mass increase in the height range from h; to hp. 


leave the COR FOV, the measured mass is lower than this 
linear-model prediction. 


4. SOURCES OF CME APPARENT MASS GAIN 


As mentioned before, the CME apparent mass might consist 
a significant contribution from ambient compressed solar wind 
plasma. In some appropriate circumstances, shock might form 
around a CME moving faster than the fast magnetoacoustic 
wave. For events E3, E5, in COR images we saw deflected 
streamers which were not adjacent to the CME, and left a very 
diffusive space between the deflected streamer and the CME. It 
might be the signature of a shock (Vourlidas et al. 2003; 
Ontiveros & Vourlidas 2009). It is almost impossible to 
separate the brightness contribution of the shock from that of 
the CME. Although solar wind pileup and compression due to a 
MHD shock are different physical processes, as the observed 
shocks were far fainter than the CME, we did not remove these 
possible shocks from our mass analyses and treated them as an 
extreme case of compression. 

To quantify the piled-up mass of the solar wind plasma 
surrounding the CME, a highly ideal model, called the “snow 
plow” model, is applied. In this model, a CME accretes mass as 
it propagates and its momentum is changed through the 
interaction with the ambient solar wind. The mass change per 
unit time is described by (see also in Tappin 2006) 


dm 
— = pA |Vsw — Vv 


dt i 6) 


where m, v and A are the mass, speed and the cross-sectional 
area of a CME, respectively, and vw and p,w are the ambient 
solar wind speed and density, respectively. The value of A is 
approximated by the solid angle Q in Table 1 times h’, i.e., 
A = 46(a + 6)h?. The derivation of the above equation could 
be found in the Appendix. It should be noted that the equation 
assumes that the entire piled-up plasma moves together with 
the CME. Although this assumption may not be true, it is still 
worth investigating if this model can represent the mass gain 
during the CME propagation. Actually the “snow plow” 
equation gives the maximal amount of mass that the solar 
wind can contribute to CME mass. 


4.1. Parameters of Background Solar Wind 


The “snow plow” model requires the solar wind parameters. 
Thus we use 3D MHD simulations to obtain the undisturbed 


ms 


zm; Measures the contribution of the solar wind pileup to 
m2? — m 


solar wind parameters along the CME path from 2.5 to 20 Rs. 
The solar wind density and speed at different heights are 
calculated as the average over a spherical sector of 60° which is 
centered at the 3D propagation direction of a CME. The 
numerical scheme in the 3D MHD simulations is a total 
variation diminishing (TVD) scheme in a Sun-centered 
spherical system (Feng et al. 2003, 2005; Shen 
et al. 2007, 2009). For all six events, the synoptic maps of 
longitudinal magnetic field over a Carrington rotation from 
Wilcox Solar Observatory are used as input to the code for 
reality. The details of how to get a steady background solar 
wind could be found in our previous works (e.g., Shen 
et al. 2013b; Wang et al. 2014). 

The averaged solar wind density and speed against height 
over the spherical sectors are presented in Figures 8 and 9, 
respectively. For all six events, the variation of density follows 
about h © to h-* below 7-8 Rs, and follows h`? above 7-8 Rs. 
It is also the height where the occulter effect starts to be 
negligible. The solar wind speeds have a rapid increase at lower 
heights and then a slower increase at larger heights. The speeds 
have a similar trend of evolution for different Carrington 
rotations, but reach different values at 20 Rs which range from 
about 220 to 460kms~'. The solar wind parameters at 20 Rs 
have been listed in Table 2 for reference. We also plot in each 
panel the speed of the CME derived from a quadratic fit to the 
height-time diagram. The scale of the solar wind speed is 
marked on the left Y axis, and the scale of the CME speed is 
marked on the right Y axis. For all events, we find that 
Veme — Ysw decreases with height. 


4.2. How Much Mass Can Solar Wind Pileup Contribute? 


The calculation of the mass that the solar wind pileup can 
contribute is straightforward, as all the parameters m, v, A, Psw 
and vw can be obtained from either CME observations or 
MHD numerical simulations. The CME speed v was estimated 
from a quadratic fit to the height—time plot, as it provides a 
better result than the linear fit. We started the mass pileup from 
m; at hı, below which we assumed that the mass increase is due 
to the geometric occulting effect and the mass supply from 
lower corona. If the mass increase above h, was solely because 
of the solar wind pileup, we would obtain the mass evolution 
indicated by the red line in each panel of Figure 7. The shadow 
region in orange is obtained by assuming that the solar wind 
density and velocity have an uncertainty of 20%. 

It is found that, for most events (E3 through E6), the solar 
wind pileup mass is not sufficient to interpret the measured 
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Figure 8. Solar wind number density against height derived from MHD simulations. The number density is the average of the density over a cone of 60° centered at 
the CME propagation direction. Three dashed lines indicate the variation of density following h”, h~*, and h~°, respectively, from upper to bottom. The vertical 


dotted lines mark the height at 7.5 Rs. 


mass increase from COR observations represented by black 
symbols, suggesting that apparent mass increase is also 
contributed by the outflow from the low corona. For the event 
El, the solar wind pileup mass has a higher value than the 
measured ones above h; in the late-phase evolution in the COR 
FOV, while for the event E2, it is more or less consistent with 


the measured mass. Note that the “snow plow” model gives an 
upper limit of the mass estimate. This model assumes that the 
solar wind plasma colliding with the CME will all be attached 
to the CME. However, this is an ideal case. Some solar wind 
plasma may not be attached, and just flow around the CME and 
eventually become the ambient background. Therefore, the 
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Figure 9. Solid lines: the solar wind radial speed against height derived from MHD simulations, which is the average of the speed over a cone of 60° centered at the 
CME propagation direction. Plus signs: the CME speed. 


actual mass increase due to the solar wind pileup may drop 


below the red lines. 


Although the solar wind pileup mass may be overestimated, 
the results still imply that the solar wind pileup probably makes 


a significant contribution to the apparent mass increase 
observed by the coronagraphs. By comparing the piled-up 
mass m, from h; to h2 with the total mass gain mz — m during 
the same period, we find that the piled-up mass occupies more 
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than half of the total mass gain as listed in the column for f, in 
Table 2, or m, is at least 19% of the final mass m, before 
the CMEs leave the COR FOV. Since the piled-up mass 
are overestimated and our infinitesimal-angular-width assump- 
tion may make the measured mass underestimated, the ratios of 
the piled-up mass to the total mass gain given above are 
obviously overestimated. However, even if we consider a 
100% overestimate, the ratio of the piled-up mass to the total 
mass gain is larger than 25% in the height range from h; to hy. 
That is to say, although the pileup is not sufficient to explain 
the observed mass increase, its contribution is non-negligible. 
Due to the overestimation, it is not clear whether the pileup 
could make a major contribution to the measured mass 
increase. 

Besides, from Figure 7, we can find that the slopes of 
mass evolution dm/dh of the measured mass and 
modeled pileup mass gradually decrease with height. When 
the CME was about to leave the COR FOV, the two slopes 
are nearly the same for most of the events except the first one. 
It indicates that solar wind pileup makes a more important 
contribution to the mass increase at larger distances from 
the Sun. 


5. CONCLUSIONS AND DISCUSSIONS 


Based on the brightness enhancement in coronagraph 
images, we have followed the mass evolution of six CMEs. 
The deprojected speed of the investigated CMEs covers from 
450 to 1300kms'. All of these CMEs with different 
kinematics have an increasing trend in their mass evolution 
even when the occulting effect is removed. The long-lasting 
accumulation of CME mass in coronagraph FOV implies that 
the initial kinetic energy in the CME source region is smaller 
than the energy estimate using the final asymptotic mass 
(Carley et al. 2012; Feng et al. 2013b). 

Physically, there may be two sources of CME mass gain in 
the corona: the solar wind mass piled up around the CMEs and 
the mass supply by the outflow from the dimming regions in 
the low corona (Aschwanden et al. 2009; Jin et al. 2009; Tian 
et al. 2012). We calculated the mass contribution of the solar 
wind pileup from the height beyond which the occulting effect 
is negligible. It is found that solar wind pileup may make a 
non-negligible contribution to the apparent mass increase 
observed by coronagraphs. For all of the events, the solar wind 
piled-up mass might occupy more than half of the total mass 
gain during the same period from 7 or 8 Rs to the edge of the 
FOV of COR2. It has to be pointed out that those ratios 
represent upper limits of the pileup contribution. However, 
even if we consider a 100% overestimation, the ratio of the 
pileup mass to the total mass gain could be larger than 
one-fourth in the height range from 7 to 15 Rs. We also find 
that the contribution from the solar wind pileup becomes 
increasingly significant as a CME propagates to larger 
distances from the Sun. In short, Our work reveals that the 
solar wind pileup is not sufficient to explain the measured mass 
increase in the COR FOV. However, its contribution is non- 
negligible. Due to the overestimation, it is not clear yet whether 
the pileup could make a major contribution to the measured 
mass increase. 

Whether the solar wind pileup is a major contributor to the 
final apparent mass in the COR FOV is not clear either. If we 
assume that the mass below 7 or 8 Rs came from the mass 
outflow from the dimming region, the ratio of m,/mz reveals 
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that the solar wind pileup comprises less than half of the final 
mass for all events. In this case, it may not be a dominant 
contributor to the CME apparent mass in the COR FOV. 
Whether this assumption is valid or not requires further 
investigation. Future work involving the mass estimate from 
the dimming regions using SDO/AIA multi-wavelength data 
will be pursued. If the solar wind pileup could also contribute 
to the mass below about 7 Rs, for some of the CMEs, the pileup 
may be a dominant contributor. 

Both the piled-up mass derived from the “snow-plow” model 
and the “virtual mass” introduced in Cargill (2004) are related 
to the solar wind mass density. However, they are theoretically 
and observationally different. From a hydrodynamic point of 
view, the “virtual mass” is the “added mass” of an accelerating 
body moving in a fluid, and is required to obtain the correct 
accelerating force in the momentum equation. The “virtual 
mass” in Cargill (2004) was derived under the assumption of an 
incompressible fluid, which means that the density in the 
surrounding solar wind does not change over time. The total 
brightness observed with a coronagraph is proportional to the 
density, as there is no density change caused by the added 
mass, we will not be able to detect it from base-difference 
images. On the other hand, the piled-up mass in the “snow- 
plow” model is a compressional effect, and the resulting 
density enhancement could in principle be detected in the base- 
difference images. Furthermore, pileup occurs even when a 
CME is moving steadily without any acceleration. 

DeForest et al. (2013) tracked the CME mass from 
STEREO/COR2 FOV to HI2 FOV, and found that the mass 
of the CME under investigation continuously increased until it 
reached a distance of about 0.7 AU, but the mass increase 
dropped with distance. It is consistent with the trend of 
evolution predicted from our mass evolution in the COR FOV. 
When the CME started to leave the COR FOV, the measured 
mass increase per height (dm/dh) seems to agree with the 
prediction of the “snow plow” model. As an outlook of the 
work in this paper, we will test whether this model can explain 
the mass evolution in the HI FOV. Furthermore, if we assume a 
proper shape of CME at 1 AU, the mass residing in the CME 
and in the shock sheath can be approximately calculated from 
in situ measurements, which can add another data point in the 
mass evolution profile. 
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APPENDIX 
“SNOW PLOW” MODEL FOR CME PROPAGATION 


“Snow plow” model assumes that all the solar wind plasma 
colliding with the CME will be attached to the CME. Based on 
this assumption and considering the 1D problem, the mass 
exchange between the CME and the solar wind is given by 


Am = —Amw = pyyA [aw — v|At (4) 
or 
dm dmMmsw 
L an = pA | vw — 
dt dt Psw |» i ©) 


where m, v and A are the mass, speed and the area of the cross- 
section of the CME, respectively, and my, Psw and vsw are the 
mass, density and speed of the solar wind, respectively. The 
cross-section area A = 46(a + 6)h*, where 26 is the flux rope 
edge-on angular width, 2(a + 6) is the face-on angular width, 
and A is the height of the leading edge. Using absolute value of 
(Ysw — v) means the solar wind pileup could occur at either the 
frontside or backside of the CME. When the CME speed is 
higher than the solar wind speed, the solar wind piles up mainly 
at the frontside of the CME; When the CME speed is lower 
than the solar wind speed, the pileup mainly occurs at the 
backside. For the CME events in this paper, we have CME 
speed always higher than the solar wind speed. Therefore the 
pileup mainly occurred at the frontside. 

The mass exchange causes the loss of the momentum of the 
background solar wind, which is 


AP,w = AlswVsw = — PywA | Reve _ v| vewAt (6) 
or 
dp dm, 
E ~ a Vsw = —PswA [vsw = v[vsw: 0) 
Meanwhile, the conservation of momentum requires 
dp APs i 
z ee 
dmv dmsw 
——— = — - Vy, 8 
dt dt 2 
dm, 
>—v + ma = = sw (9) 
dt 


where a is the acceleration of the CME. Using Equation (5), the 


above equation leads to 
dm 
ma = — (vw — v). 


dt 
On the other hand, combine Equations (7) and (9), we have 


(10) 


(1) 


dm 
a + ma = px, AVsw | Vsw — v| 


After Substituting Equation (10) for ma in the above equation, 
it is obtained that 


dm 


12 
i (12) 


= Pw å | Pee — v| 
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which is the equation we used to calculate the mass 
contribution from solar wind pileup. 
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